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From a simple Lagrangian the equations of motion for a particle with spin are derived. 
The spin is shown to be conserved on the particle's world-line. In the absence of a spin 
the equation coincides with that of a geodesic. The equations of motion are valid for 
massless particles as well, since mass docs not enter the equations explicitly. 

1. Introduction 

According to the equivalence principle (EP) of the general theory of relativity 
(GR), motion of structureless test particles in a gravitational background is de- 
termined only by the spacetime geometry: particle worldlines are the geodesies of 
the spacetime. Things become more complicated for test particles which are not 
structureless and carry, for example, a non-vanishing charge or spin. In such cases, 
the worldline of a particle, in general, is no longer a geodesic, but is modified by 
electromagnetic and/or spin-gravity forces (for, e.g., [1, 2] and references therein). 
The problem of the motion of classical spinning particle in external fields has oc- 
cupied scientists for nearly all of the last century. Covariant equations of motion 
for a relativistic particle in an electromagnetic field were first written more than 
70 years ago by Frenkel For the case of a gravitational field, Mathisson * and 
Papapetrou ^ found the energy-momentum and angular momentum propagation 
equations for a rotating test body ( "pole-dipole particle") according to the Ein- 
stein's GR. Tulczyjew ^, Beiglboeck '' and Madore ^ developed these into laws of 
motion by adding a definition of a centre-of-mass world line. Later, Dixon ^ gener- 
alized these treatments and made them more rigorous. All theses works dealt with 
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the motion of extended rigid bodies or tops. It is far from obvious whether one 
can observe in practice the spin corrections to the equations of motion of elemen- 
tary particles. However, the problem of influence of the spin on the trajectory of a 
particle in an external field is not of only theoretical interest. Spin-dependent cor- 
rections certainly exist in differential cross sections of scattering processes. It was 
proposed long ago to separate charged particles of different polarizations through 
the spin interaction with external fields in a storage ring of accelerator Though 
this proposal is being discussed rather actively now (see a review in Ref. [11]) it is 
not yet clear whether it is feasible technically. The EP can be put to test in an as- 
trophysical setting, a recent proposal being based on the analysis of the differential 
time delay between the arrival of left and right-handed circularly polarized (LCP 
and RCP) signals from the millisecond pulsar PSR 1937+214 However, far few 
papers exist on the theoretical foundations of possible deviations of the photon mo- 
tion from geodesies (see, e.g., Ref. [12, 13]). Here we report yet another approach 
to this long-standing problem. If an appropriate Lagrangian density is taken into 
account, then the photon equation of motion (modified geodesic equation) can be 
found from the Euler-Lagrange equations. 

2. Formalism 

The generalized concept exists that classical particles follow the path of the least 

spacetime distance between the endpoints, even when space is curved by gravity. 
Thus, in (pseudo-) Riemannian spaces the geodesic equation is found from the 
variation of an action S, identified with the parameter s of a curve interpreted as 
its length. The same method has been used to find the geodesic equation of light 
where, however, one technical problem arises: photon's worldline has no length. One 
way to avoid this difficulty is to consider the motion of a massive vector particle 
and, if the obtained equation does not contain mass explicitly, simply put the mass 
to zero; with one further condition that the four velocity of a particle be a null 
vector. This method has been widely used for scalar particles and is known to give 
the equation of geodesic regardless of mass. In this work we apply this approach, 
which makes it possible to use the length parameter of the action principle. 

Usually in order to describe the behavior of a field in a given gravitational 
background, one solves the corresponding field equations for a given metric. If 
the goal is to describe waves, one can take the corresponding wave solution. In a 
spherically symmetric spacetime the solutions contain factors expressed in terms of 
spherical harmonics. This method works well when the wavelength is comparable 
to the scales under consideration. It is not so in the case of light propagating in 
the vicinity of a massive object and, to describe the propagation of light as an 
electromagnetic field, one would have to employ the spherical harmonics of a very 
high order. The corresponding solution would look too complicated and tell little 
about the behavior of light. 

Many efforts have been spent in the last decades to work out a simple approach 
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to this problem that would give a satisfactory approximation to the wave as some 
curves that could be called "rays" and, at the same time, would take into account the 
polarization of light. In the present work we try to work out a simple approximation 
of this type. Our idea consists of the following: we consider a massive vector 
field obeying the Proca equation, describing the propagation of this field, in some 
restricted domain of the spacetime. The shape of the domain can be chosen as that 
of some world-tube transverse to the wave, with the cross-section comparable to 
the wavelength, or, say, not greater than two orders of magnitude. As this tube is 
timelike there must exist a timclike curve A in its interior, which specifies a local 
coordinate system as the time axis. And if the tube is not too wide, this coordinate 
system would cover the entire interior. If s is the proper time on the curve A, the 
curve can be chosen in such a way that the field equation reduces to 

D DAi , , 

m^A, = Q, (1) 



Ds Ds 



the same way as it happens in standard Cartesian coordinates for Minkowskian 
spacetime. Correspondingly, the field Lagrangian is — A^ + m^A^, where dot stands 
for the covariant derivative on s, if and only if the curve A is chosen properly. This 
Lagrangian should contain one more term responsible for the shape of A, yielding 
the geodesic equation when we switch off the field. The form of this term is well 
known: l/2mx^, thus our final Lagrangian is 

2L = mx^ - A^ + m^A^ , (2) 

and coupling between the field and the shape of the curve A is incorporated in the 
form of a covariant derivative A, which contains the product of connection, velocity 
X and the field. 

The derivation of the conservation laws is more convenient in orthonormal 
frames. In what follows, e" will denote the components of an orthonormal 1-form 
frame field, 

= e1{x)dx' , (3) 
and ejj the components of its dual vector frame field, 

ea = 4(a;)^- (4) 

Here frame indices are always a, b, c, coordinate indices are i, j, k, .... The metric 
tensor can be expressed as 

g = g^jdx" (g) dx' = SabO" «> 9^ . (5) 

The connection 1-form for these frames may be introduced through the first struc- 
ture equation: 

dO" = 1^^" A6\ UJab+^ba=Q (6) 
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and the connection coefRcients •jabc are that of the expansion of this 1-form in the 
local frames {9°'}: 

LO,^=lab0''. (7) 

Thus, 

or A"' is a covariant derivative in orthonormal frame on a curve x^{s) with 7^^" 
being a spin connection. 

3. Equations of Motion 

Each generalized coordinate has its conjugate generalized momentum: 

dL 

P<^ = ^= 'mVabi'' - r]dbA''A''jJ . (9) 



and 
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Corresponding to them the Euler-Lagrange equations are 
and 

Let us first consider equations for the Ea (12). Right hand side of this equation is: 

1^ = ~A^ri,,^Jx^ + m'A^rjae ■ (13) 
And Euler-Lagrange equations for Ea will be: 

^Ea - EaiJx^^ - m^A'^rjac = , (14) 

where we used Eq. 10. Expression in the brackets is a covariant derivative for Ea, 
thus, we obtain: 

Ea - m'^A'^rjac = . (15) 
Again using Eq. 10, we finally obtain: 

A%ab + m^A^ab = , (16) 
which, as we can see, reduces to Proca equation for the four- vector field Ai^{x) ^^: 

(□ +m^)Ap = 0. (17) 
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We describe spin of the particle in this model directly by a tensor of spin Sab- 

1 



mSab = 2 [AhAa - AaAb 

To obtain the equation of motion for spin, we write 

dSab ^ ^ [ A A i A 

'^~ds' 2 (is V ~ * " 
Due to (17) all partial derivatives of A vanish and we are left with: 



(18) 



(19) 



dS, 



m- 



db 



ds 



A{Aj AjAi 



d_ 
dx 

AaAb AbAa 



AdAa AaAd 



(20) 



Here is the matrix introduced in the equation (4) and derivatives (e^) are 
obtained from the first structure equation (6). 

Using our definition of a spin tensor (Eq. 18), we obtain: 



dS, 



db 



ds 



or 



X^ilcdSab + lcb'Sad) , 

DSdb 



Ds 



= . 



(21) 
(22) 



Thus, spin is transported parallel to itself along the worldline. 

To derive Euler-Lagrange equations for generalized momentum, we return to the 
coordinate basis. In this case Lagrangian (Eq. 2) will take a form? 



dA^ 



2L = mgijx'x' - gij { ^ + Tlix'^A^ 



fdA> 

\dr 



TiiX^A'\+m'gijA'A\ (23) 



where Fj.^ now are Christoffel symbols. Generalized momentum from here is: 

dL 



In coordinate basis Eq. 11 becomes: 

d 



ds 



Pi 



dL 
dx^ 



LHS of this equation is: 



S = ^ + mx^i'^dkg^j - x^A'AmdjT^ - 



ds ds 



(24) 



(25) 



, (26) 



"It must be noted that while Eq. 2 represents a non-relativistic Lagrangian of the free mass point 
when we remove the field part, here we use the relativistic form; in this case both forms are 
actually equivalent. When we are considering relativistic kinematics (free motion of a mass point, 
for example), the trajectory will be a geodesic regardless of what form a Lagrangian we take. 
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where derivatives di are defined as di = RHS is: 

^ = \mdigjkx'x'' + diQmn ("i^i" + m^A-A") - gmnx" A' d^T^i . (27) 
Subtracting (27) from (26), we obtain: 

dx^ 1 dA"^ 

^9ij-^ + mx^x''dkgij - -mdiQjkX^x'' - digmnni^ A"^ A"" - r^A'-^ + 

d 

digr^nA^^A^ - ^uAm^ - x^A'A^d^T^ + g^ni'' A' A^ diVf, = . (28) 

Using 

drgmn = \{gknT''^i+gkmT^in) , (29) 

and some algebra, we obtain from (28) 

(2k I " \ k ' I ' ' ' 

"^^'"'"Ds v"^ ^miA"^Ak + T'j^A Am j — V^T mjX-' A Ak + V^^A™" An — 

^ii A Am + x'^ A Am^ii — x'' A An^djTj^i -\- x A AmdiTi^i = . (30) 

We may notice that terms in the brackets are zero due to the Proca equation (17) 
and terms containing AA cancel; the remaining terms 

gi,w.^ = x^A^A, {djT^, - diV'^, + rtrj} - r^,r^) . (3i) 

With the definitions of the curvature tensor and spin tensor (18), the equation takes 
the evidently covariant form: 

9ij^=R'' jui'S,^ . (32) 

This equation coincides with the Papapetrou equation ^. It must be pointed out 
that in case of a zero spin this equation becomes a geodesic. If we reparametrise 
the curve with some new parameter A in such a way that 

dx^ dxP 



we can rewrite Eq. 32 as 

Dx^ 

"da 

This equation is valid for the massless particles with spin as well, since mass does 
not enter the equation explicitly. 



9ii^^=R'' jiii'S^' . (34) 
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